The Neutron Peak in the Interlayer Tunneling Model of High Temperature 

Superconductors 
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Recent neutron scattering experiments in YBCO exhibit an unusual magnetic peak that appears 
only below the superconducting transition temperature. The experimental observations are ex- 
plained within the context of the interlayer tunneling theory of high temperature superconductors. 



PACS: 74.72.-h, 74.20-z, 74.25Ha 

Many experimental observations in high temperature 
superconductors are commonly fit with a phenomeno- 
logical model that derives from the original theory of 
Bardeen, Cooper and Schrieffer (BCS). The model has 
certain ingredients. It is characterized by a gap equation 
corresponding to a presumed symmetry of the order pa- 
rameter, with an adjustable dimensionless coupling con- 
stant, or alternately an adjustable ratio of the gap to 
the superconducting transition temperature, and a given 
Fermi surface. Inherent in this description are the co- 
herence factors that determine a number of interference 
processes unique to the BCS theory. Although a mi- 
croscopic derivation of this effective model for the high 
temperature superconductors does not exist, the model 
is still used with a considerable degree of confidence. The 
main difficulties, to which we return below, are the un- 
usual normal state properties of these materials and the 
enormously high transition temperatures, but there are 
many others. Such difficulties are extensively surveyed 
in the literature |l]]. 

In the absence of a microscopic derivation, it is useful 
to ask if this phenomenological BCS model is unique and 
if an alternative phenomenological model exists that is 
capable of capturing features of these superconductors. 
One such physically motivated model was elaborated in 
a recent paper j|. In the present paper we examine this 
model, called the interlayer tunneling model, to interpret 
the startling neutron scattering experiments in optimally 
doped YBCO @. 

These neutron scattering experiments are startling for 
a number of reasons. The experiments appear to have 
established that there are no sharp, or even broad, fea- 
tures in the magnetic excitation spectrum in the normal 
state. In contrast, the superconducting state exhibits a 
sharp magnetic feature that appears to be localized both 
in energy and momentum. It is located at an energy of 
41 meV and near a wavevector (n/a, n/a, w/cb), where a 
is the lattice spacing of the square-planar CuO lattice, 
and Cb is the distance between the nearest-neighbor cop- 
per atoms within a bilayer. The intensity under the peak 
vanishes at the transition temperature, but its frequency 
softens very little. 

In principle, there can be many explanations, and such 



explanations have been proposed [Q . None of these ex- 
planations are fully microscopic, nor are they fully con- 
sistent with some of the other prominent experimental 
observations in these materials. In particular, we draw 
attention to the fact that explanations that rely on form- 
ing coherent linear superpositions of the bilayer bands 
would also predict large splittings of the bands in certain 
regions of the Brillouin zone. Such splittings, studied 
extensively in experiments, are missing to a high degree 
of accuracy ||. The argument that the observed bands 
are highly renormalized has little force. If the renormal- 
ization effects are so strong as to reduce the splitting 
to nearly zero, it is not meaningful to speak of a coher- 
ent superposition of states at the experimentally relevant 
temperatures. 

The explanation we have to offer within the interlayer 
tunneling model is exceedingly simple and purely kine- 
matic in origin. We also present illustrative but quanti- 
tative results to compare with experiments. The compar- 
ison shows that many of the features of the experiments 
are captured well. While many improvements can possi- 
bly be made, the kinematic aspects should survive future 
scrutiny. 

The model Hamiltonian, motivated earlier is that 
of a bilayer complex. Its generalization to many layers is 
straightforward and not essential for the present purpose. 
It is 
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Here i = 1,2 is the layer index. The fermion operators 
are labeled by the spin a and the wavevector k; V is the 
in-plane pairing interaction. The last term describes tun- 
neling of pairs between the layers. Such a Hamiltonian 
should be understood from the point of view of an effec- 
tive Hamiltonian, as emphasized recently p] . It incorpo- 
rates the unique feature of the interlayer mechanism, that 
tunneling occurs with conservation of transverse momen- 
tun k. Therefore, the momentum sum in the Tj term is 
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only over k and q. Disorder between the layers is weak, 
and even disorder would not affect this crucial difference 
between the tunneling and interaction terms. 

Only in the subspace in which both the states (k f) and 
(— k |) are both simultaneously occupied or unoccupied, 
is the following reduced Hamiltonian sufficient: 
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This is because these are the only matrix elements that 
survive. It must be remembered that the actual Hamil- 
tonian is H and that -ff re d is merely a convenient tool 
to calculate the matrix elements in the paired subspace. 
We do not anticipate that the Schrieffer zero-momentum 
pairing hypothesis will be necessarily as accurate in this 
case as in BCS, but it surely will be a guide. 

An important feature is the absence of the single par- 
ticle tunneling term from the Hamiltonian, which allows 
for a coherent superposition of the states of the layers 
that leads to a splitting, not observed in experiments. 
Therefore, such terms are not included in our model. 
Although the single particle tunneling term is absent, 
incoherence does not preclude a particle-hole two parti- 
cle tunneling term |^|| , because the two particle term is 
generated by a second order virtual process and the ques- 
tion about coherence is irrelevant; the energy need not be 
conserved in a virtual process. However, this term does 
not appear to be crucially important for optimally doped 
materials (except, see later). For underdoped materials 
exhibiting "spin gap" phenomena (?j , this neglected term 
can be quite important ||. We also note that we shall 
treat the Fermion operators as ordinary anticommuting 
fcrmion operators. This can only be a crude approxima- 
tion due to the non-Fermi liquid feature of the normal 
state. A possible improvement along the lines discussed 
earlier is exceedingly complex. 

The mean field treatment of the reduced Hamiltonian 
is straightforward B and leads to the gap equation: 
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where Xk = (Ak/2i7k) tanh(i?k/2T) is the pair suscep- 
tibility, with E k = y/(e k - fi) 2 + A k . Here /i is the 
chemical potential. Note that, until now, we have not 
specified the symmetry of the in-plane pairing kernel. 
We shall now assume it to be of the symmetry d x i_ y i 1 
Vk,k' = Vg k g k >, where g k = \ [cos(fc x a) - cos^a)]. 

Further specifications are necessary. The in- 
plane one electron dispersion is chosen to be e k = 
—2t [cos(k x a) + cos(k y a)] + it 1 cos(k x a) cos(k y a). The 
quantitative calculations are carried out with a repre- 
sentative set of parameters appropriate to YBCO. These 



are t = 0.25 eV, t' = 0.45i, and /z = -0.315 eV, cor- 
responding to an open Fermi surface, with a band fill- 
ing of 0.86. The quantity Xj(k) was argued to be 
T 7 (k) = (Tj/16) [cos{k x a) - cos(fc y a)] 4 . The validity of 
this expression is now strengthened by detailed electronic 
structure calculations || for the single particle hopping 
matrix element ij_(k). 

The magnetic neutron scattering intensity is propor- 
tional to the imaginary part of the wavevector and the 
frequency dependent spin susceptibility x(q, u>), which 
for the above model is simply the expression: 
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n k>q H = w - (£ k+q - E k ) + is, n£± (w) = u ± (£ k+q + 

E k ) + iS, F± q = /(£ k +q)±/(£k), and f{x) is the Fermi 
function. 

Note that near T — only the A~ terms contribute, 
and these are negiligible unless Ak and Ak+ q are of op- 
posite sign, as noted in Ref. In fact, nothing is ob- 
served for q w or 2ir, as expected, and the peak appears 
at q = (tt/cl, tt/o), connecting two opposite-sign lobes of 
x 2 — y 2 . At higher temperatures, the A + part is tem- 
perature dependent, but almost frequency independent 
function for experimentally relevant frequencies. The \ 
with the A + term omitted will be denoted by x- We shall 
first calculate the imaginary part of this susceptibility by 
restricting all wavevectors to be in the CuO-plane. Later, 
we shall discuss the dependence on the momentum trans- 
fer perpendicular to the plane. 

At T = 0, Im x is calculated by solving the gap equa- 
tion for a set of Tj, and with A = N(0)V = 0.184, where 
iV(0) is the density of states per spin at the Fermi energy 
(V = 0.2 eV, N(0) = 0.92). These are shown in Fig. 
While the calculations show a peak at the threshold for 
Tj ^ 0, in the pure BCS case, Tj = 0, the intensity 
is a step discontinuity at the threshold, and other exci- 
tonic enhancement mechanisms are necessary to produce 
a peak at the threshold. 

For illustrative purposes the gap equation is solved at 
finite temperatures for A = 0.184 and for a fixed Tj = 
0.075 eV. No attempts were made to fit precisely the 
experimental data. The results are plotted in Fig. ||. 
We calculate the intensity under the peak and plot it as a 
function of temperature. The intensity was calculated by 
fitting a linear background at high energies. The results 
are shown in Fig. ||. The intensity falls off to zero as the 
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FIG. 1. Im x at T = 0. From left to right, Tj =0.025, 0.05, 
0.075 eV. The curve corresponding to a step discontinuity at 
the edge is for BCS, with the maximum d-wave gap equal to 
0.025 eV. 
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FIG. 2. Im x for Tj = 0.075 eV. From left to right, T =120, 
110, 100, 80, 60, 40, 20, and K. The results for T = and 
20 K are almost indistinguishable. 



temperature is raised to T c and the results compare well 
with experiments. Also shown in Fig. [?], is the position 
of the peak, whose softening is weaker than the fall-off 
of the intensity. However, the experimental softening of 
the position of the peak is even weaker. 

The question to address is why Im x produces a peak 
in the interlayer tunneling model, while the BCS model 
exhibits only a step discontinuity at the threshold. The 
reason is simple kinematics that holds for an open Fermi 
surface combined with an unusual feature of the inter- 
layer gap equation. This feature was emphasized in early 
works on the gap equation 10 1, but has not been specif- 
ically noted in recent calculations: it is that when Tj 
dominates, as is necessarily so for high T c , the density 
of states for a given value of Tj/2 is approximately a 6- 
function at Tj/2. Concomitant with this is the fact that 
electrons with \sk — (i\ > Xj(k)/2 are unaffected by pair- 
ing, in contrast to BCS where effects of the gap extend 
out to high energies. 

For simplicity, consider T = 0, set the coherence fac- 
tor to unity and look at the wavevector Q = (it /a, n/a). 
In the BCS case, that is with Tj(k) set equal to zero, 
(£k + ^Q+k) has cither minima, or saddle points. The 



FIG. 3. The intensity (solid triangles) and the position of 
the peak (filled circles) normalized to the zero temperature 
values. The position of the peak is at 41.2 meV at T = 0. 



minima give rise to a step discontinuity at the threshold 
at approximately twice the maximum of the <i-wave gap 
(See Fig. [j].). The saddle points give rise to van Hove sin- 
gularities which are generally not at the threshold, but 
can be brought close to the threshold by adjusting pa- 
rameters. 

Now consider the case of very small in-plane pairing 
kernel in the interlayer tunneling equation. The unusual 
feature is that the superconducting region is a very nar- 
row region around the Fermi line as shown in Fig. |]. The 
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FIG. 4. The superconducting region and its mapping under 
k — > Q — k. The diamond shaped overlap regions contribute 
to the peak of the imaginary part of the spin susceptibil- 
ity. The dashed arcs are the constant aj-contours centered at 

(M). 

scattering surface is now radically different. The contri- 
bution to the peak must come from regions in which both 
Ak and Ak+Q are finite. Also, note that £k+Q = -^Q-k- 
The image of the superconducting region including the 
Fermi line under the mapping k — * Q — k is shown in 
Fig. [§ The regions in which both Ak and Aq_k are fi- 
nite are the diamond shaped overlap regions, with one of 
its diagonals given by k x + k y = n. Consider the lower di- 
amond. In this region, we can write i?k ~ Tj cos 4 (d/\/2), 
where d is the distance from (ir/a, 0). The corresponding 



3 



constant w-contours are the arcs shown in Fig. [|. 

Note that the Im \, as a function of u, is proportional 
to the length of the arcs for which (E^ + Eq^) = w - The 
lower and the upper cutoffs are defined by the arcs touch- 
ing the diamond, and the peak is contained between these 
two cutoffs. As u> increases, the arc length increases ap- 
proximately linearly within the diamond, and then drops 
approximately linearly. Moreover, because the supercon- 
ducting region around the Fermi line is very narrow, the 
width of this peak is also very narrow, as in experiments. 
The peak is at Tj(k ), where k is located at the center of 
the diamond. This elementary analysis substantiates the 
numerical results in which none of the above simplifying 
approximations were made. 

To understand the dependence of the scattering in- 
tensity on the momentum transfer perpendicular to the 
plane, it is necessary to consider the mixing of the 
electronic wavefunctions between the layers. This can- 
not be a coherent superposition of states, otherwise the 
bands will be split, which, as stated earlier, is not ob- 
served in experiments. However, a virtual mixing of 
the states due to second order processes, similar to su- 
perexchange, is possible. This virtual mixing can be de- 
scribed by approximately constructing the operators 
(ck, CT ,i + r?(k)ck, -,2)/V^ _ 

f?(k)ck,fT,i)/\/r 
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ry 2 (k), where rj(k) is the mix- 



ing parameter. The order parameter that takes the mix- 
ing into account is Ak = (a k j iO^_ k ^ 1 + a k ^ 2 q^ 



It is also possible to rewrite this in terms of the operators 
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where we have used the notation c c<! = c c ^ 
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This virtual mixing of the wave functions on the 
layers is a necessary concomitant of the source of 
the pairing energy in the hopping terms in the origi- 
nal Hamiltonian -Ek,(r,^j t -i-( k )(v,i Ck #.i + h.c.) = 
— ^ k ij_(k)(n k a — n k CT ), where n are the number op- 
erators. The pairing energy must be numerically aug- 
mented by its expectation value. However, the two- 
particle tunneling term alone does not give quite enough 
difference in occupancy to explain the experimental mod- 
ulation with q z . In a forthcoming paper JsJ , we will show 
that this tendency is strongly enhanced by including the 
"superexchange" type terms caused by virtual particle- 
hole hopping. In fact, it may be that the pairing occurs 
primarily in the "even" linear combination alone |0]. 

Then, any electron excited from the condensate must 
leave behind a hole in the even state. Similarly, it can ac- 
cept only an electron in the odd state, hence the striking 
even <-> odd selection rule observed in the neutron scat- 
tering experiments. Note that although the odd states 



are nominally unpaired, the addition or removal of an 
electron from the odd state destroys the pairing criterion 
and costs one unit of pairing energy. Thus, there is a 
gap for all one-electron excitations, even though only the 
even state is paired. 

In conclusion, we find that these neutron results give us 
a surprisingly direct and complete picture of the nature 
of the pairing responsible for superconductivity, not only 
showing us the symmetry and approximate form of the 
energy gap, but even fixing the nature of the gap equation 
and the source of the pairing energy. In particular, the 
strong correlation between the layers excludes any purely 
intralayer mechanism for superconductivity, under our 
assumption — the singlet, and very plausible on other 
grounds — that the quasiparticle pairs do not have strong 
final state interactions. 
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cussions. 



[1] 
[2] 

[3] 

[4] 



[5] 

[6] 

[7] 

[8] 
[9] 

[10] 
[11] 



P. W. Anderson, Science 256, 1526 (1992). 

S. Chakravarty, A. Sudb0, P. W. Anderson, and S. 

Strong, Science 261, 337 (1993). 

H. F. Fong et al., Phys. Rev. Lett. 75, 321 (1995) and H. 
F. Fong et al., to be published. 

I. I. Mazin and V. Yakovenko, Phys. Rev. Lett. 75, 4134 

(1995) ; E. Delmer and S. -C. Zhang, Phys. Rev. Lett. 
75, 4126 (1995); D. Z. Liu, Y. Zha, and K. Levin, Phys. 
Rev. Lett. 75, 4130 (1995); N. Bulut and D. J. Scalapino, 
Phys. Rev. B 53, 5149 (1996); B. Normand, H. Kohno 
and H. Fukuyama, J. Phys. Jpn. 64, 3903 (1995). 

H. Ding et al, Phys. Rev. Lett. 76, 1533 (1996); D. Van 
der Marel et al. Physica C 235-240, 1145 (1994) and 
references therein. 

L. Yin and S. Chakravarty, Int. J. Mod. Phys. 10, 805 

(1996) . 

A. G. Loeser et al., Science, in press, H. Ding et al., 

Nature, in press, and references therein. 

P.W. Anderson, preprint, in preparation. 

A. I. Liechtenstein, O. Gunnarsson, O. K. Andersen, and 



R. M. Martin, preprint ( |:ond-mat/950910l| ). 
P. W. Anderson, Physica C 185-189, 11 (1991); J. Phys. 
Chem. Solids 54, 1073 (1993); Chapter VII of a forth- 
coming book. 

If we augment our Hamiltonian (Eq. |^) with an equal 
particle-hole pair hopping (superexchange) term and ig- 
nore the small in-plane interaction, it acquires an addi- 
tional approximate symmetry with respect to the inter- 
change of particles of a given spin species and of wave 
vector k between the layers. The odd order parameter 
does not respect this symmetry and must vanish unless 
the symmetry is spontaneously broken. 
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